Unconventional Mott Transition in K a; Fe 2 _,,Se2 
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Whether the newly discovered K x Fe2-i/Se2 systems are doped Mott or band insulators is key to 
how superconductivity emerges at lower temperature. With extant theoretical studies supporting 
conflicting scenarios, a more realistic approach is urgently called for. Here, we use LDA+DMFT to 
study this issue in K x Fe2- !/ Se2. We find that the undoped KFei.eSea system is a new kind of Mott- 
Kondo insulator (MKI). Electron doping this MKI drives a Mott transition to an orbital-selective 
non-Fermi liquid metal. Good agreement with spectral and transport responses supports our view, 
implying that superconductivity arises from a doped Mott insulator, as in the high-T c cuprates. 



Recent finding of high-T c supercxonductivity (HTSC) 



with T c > 30 K in K x Fe 2 - y Se 2 || and (K^Fe^Ses @,[| 
is remarkable for the following reasons: (i) While its 
SC T c is comparable to that for the 1111- and 122-Fe 
pnictides (FePn), it occurs near antiferromagnet (AF) 
insulators: Fe composition is the tuning parameter for 
insulator-metal and SC instabilities. Evidence for Fe va- 
cancy order Q is interesting: its generic effect is to re- 
duce LDA (local-density-approximation) bandwidths @ . 
Theoretically, both Mott [1-13] and band [§] insulating 
states have been proposed as candidates, (ii) It inten- 
sifies the fundamental debate [9( on the degree of elec- 
tronic correlations in Fe-based SC. Transport and optical 
data reveal insulating behavior (albeit with small gap) 
well above Tn(S), S being the electron doping. However, 



optical [lOj and ARPES studies also clearly show 



large-scale spectral weight transfer (SWT) (over an en- 
ergy scale 0(2.0) eV) as a function of temperature (T) 
across the magnetic and SC instabilities, a fingerprint of 
Mottncss. Both p(T) ~ T and bad metallicity above T c 
are features shared along with other non-Landau-Fermi- 
liquid (LFL) metals close to a Mott or selective-Mott 
instability [12]. ARPES 11] also seems to show a curi- 
ous co-existence of Mott- and "band-insulating" spectral 
features in K^Fe2- y Se2. Large-scale SWT on a scale of 
0(2.0) eV in response to a small {k B T ~O(10) meV) 
perturbation, however betrays the "hidden" strong cor- 
relations. Its implications for SC are intriguing: are the 
Mott- or "band" -like subsets of the renormalised spectra 
relevant for SC that emerges at lower T? (iii) Finding of 
large local moment value on Fe, Mf e ~ 3.3/is, suggests 
strong electronic correlations. Even in the more itiner- 
ant FePn, relevance of a "dual" picture for the parent 
magnets is now increasingly recognised [13]: such a dual 
picture must be even more relevant for K.;;Fe2-ySe2. 

These findings necessitate incorporation of reasonably 
strong multi-band electronic correlations. All five ri- 
bands crossing the Fermi energy (Ep) must be kept at 
a "minimally realistic" level in order to satisfactorily re- 
solve the doped Mott-versus-band insulator issue above. 
In multi-orbital (MO) systems, sizable correlations also 
drive new physical effects: they induce orbital selective 



(OS) bad-metallic states with no LFL coherence, natu- 
rally yielding bad metallic resistivity above T c . An OS- 
Mott scenario also generically "wipes out" a subset of 
Fermi surface (FS) sheets from the renormalized, corre- 
lated band structure (Lifshitz transition). In turn, this 
can have far-reaching consequences for the symmetry of 
the SC pair function, A(k): details and presence or ab- 
sence of gap nodes crucially depends on whether or not 
A(k) intersects the renormalized FS of such a metal. Ex- 
tant ARPES data show only electron pockets in metallic 
K0.8Fe1.7Se2 and anisotropic s-wave SC gap [3]. How- 
ever, NMR rp x (T) ~ T 2 for T < \T C remains puzzling 
in this context [15|. This unsettled state of affairs calls 
for detailed theoretical scrutiny: it must, in view of the 
fundamentally conflicting views discussed above, base it- 
self on an satisfying description of the "normal" state. 

Here, we use state-of-the-art LDA-plus-dynamical 
mean-field theory (LDA+DMFT) [lj| to address these 
issues. In LDA+DMFT studies of FePn, their degree of 
correlatedness has been the bone of contention: consen- 
sus has fluctuated between weakly [TtJ to sizably corre- 
lated limits tiaj. In K x Fe 2 _ y Se 2 , a perusal of the re- 



sistivity data [11| for non-AF but SC samples reveals 
that these are proximate to a Mott transition: Pdc{T) 
immediately above T c is very bad metallic and quickly 
crosses over to an insulator-like dependence. In these 
cases, AF order cannot be responsible for insulating be- 
havior: rather, these data show that destroying AF or- 
der reveals underlying Mottness, where (i)-(iii) above can 
be rationalised naturally. Clearly, electronic correlations 
can only get stronger as one approaches the insulator, 
and AF and/or orbital order can naturally arise via spin- 
orbital superexchange. In this work, we confirm this hy- 
pothesis, showing that a strong correlation view achieves 
good semiquantitative accord with a range of data in 
both, insulating and metallic phases of K x Fe2- y Se2. In 
particular, we clarify the co-existing Mott- and "band" 
insulator-like features in ARPES in a qualitatively new 
Mott-Hubbard scenario. Armed with these strengths, we 
qualitatively discuss the constraints our view imposes on 
mechanisms of (unconventional) SC in K x Fe2- y Se2. 
We start with the experimental structure of 
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FIG. 1: Orbital-resolved LDA density-of-states (DOS) for the 
Fe d orbitals of FeSe (dot-dashed) and KFei.6Se2 (dashed) as 
well as LDA+DMFT (with U = 4.0 eV, U' = 2.6 eV and 
Jh = 0.7 eV) for undoped KFei. 6 Se 2 . 



FIG. 2: Orbital-resolved (DMFT) imaginary parts of the self- 
energies for the Fe d orbitals of pure and electron doped 
KFei. 6 Se 2 with U = 4.0 eV, U' = 2.6 eV and J H = 0.7 eV. 
The selective-Mott nature is clear. 



KFei.6Se2 LDA calculations were performed 

using the linear muffin-tin orbitals (LMTO) [1 91 ] scheme, 
in the atomic sphere approximation: Self-consistency is 
reached by performing calculations with 242 irreducible 
k-points. The radii of the atomic spheres were chosen 
as r=2.6 (Fe), r=4.25 (K) and r=2.7 (Se) a.u. in 
order to minimize their overlap. Within LDA, the 
one-electron part is H = Ek,a, ( 7 e »W4,a, 1 T C i',» 1 <') 
where a = xz,yz,xy,3z 2 — r 2 ,x 2 — y 2 label the 
(diagonalized in orbital basis) five 3d bands. The 
MO Coulomb interactions (treated within DMFT) 
constitute the interaction term, which reads Hi nt = 



We use the MO iterated perturbation theory as an 
impurity solver for DMFT. Though not numerically 
exact, it has a proven record of recovering correct 
LFL metallic behavior [2(| and good semiquantitative 
agreement in a host of real systems. 

Fig. ED shows the LDA DOS for KFei. 6 Se 2 , whereby 
a clear and sizable reduction (O(20)%) of the average 
LDA bandwidth (Wlda) relative to that for FeSe, in- 
duced by presence of Fe vacancy order, is seen. FeSe is 
already a bad metal close to a Mott insulator [21| . and 
the significantly smaller Wlda for KFei.gSe2 then natu- 
rally suggests emergence of a Mott insulator in the latter. 
Indeed, our results show a small but clear insulating gap 
in DMFT spectra. Several interesting features, especially 
germane to the above discussion, are now manifest: (i) 
the Mott gap is clearly orbital-dependent, i.e, intrinsi- 
cally anisotropic, (ii) Examination of the orbital-resolved 
(imaginary parts) self-energies reveal a behavior hith- 
erto unique to Fe-based systems. Namely, ImE a (w) with 



a = xz, yz, xy, x 2 — y 2 clearly reveal their Mott insulating 
character, i.e, a pole at u> = Ef{= 0), as shown in Fig. [2] 
On the other hand, ImE a (w) for a = 3z 2 — r 2 simultane- 
ously shows Kondo insulator features, i.e, ImS n (w) = 
in the gap region. Thus, remarkably, KFei.6Se2 shows 
co-existing Mott and Kondo insulating gaps, and we 
dub this a Mott-Kondo insulator (MKI). One might be 
tempted to try to link this to the band insulator state 
found in LDA, especially since a Kondo insulator is an 
analytically continued version of a band insulator (where 
sizable correlations do exist above the gap scale). If 
electron doping would result in a metal where only the 
dj, z i_ r 2 band would be driven metallic, it would neces- 
sarily mandate an effective, doped Kondo insulator [22| . 
rather than doped Mott-insulator modelling. However, 
in ARPES data for the doped metal, the renormalised 
FS comprises sheets having only xz,yz orbital charac- 
ter, in accord with DMFT results (see below). Thus, the 
insulator-metal transition must now be viewed in terms 
of a doped Mott insulator, in accord with earlier model- 
based work 0, Q . We emphasise that no such insulating 
state is found for U < 3.0 eV (not shown) in absence of 
magnetic order. 

Electron doping (n = 6 + 5, with 5 > 0) the MKI leads 
to a bad-metallic state. Consistent with our results in 
Fig. [21 we find that only xz, yz orbitals show bad metal- 
lic behavior, characterised by absence of sharp LFL res- 
onances at Ep, while Mott insulating behavior persists 
in the xy, x 2 — y 2 , 3z 2 — r 2 chanels. Thus, we find an OS 
metal, and strong scattering between the Mott-localised 
and metallic states leads to complete suppression of the 
LFL quasiparticles via the Anderson orthogonality catas- 
trophe [231 ] . leading to emergence of anomalously broad 
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FIG. 3: Total DMFT DOS for KFei. 6 Se 2 (main panel) as 
a function of electron doping n = 6 + 8, clearly showing a 
low-energy pseudogap associated with bad-metallicity. Low- 
energy DMFT spectra (inset) in good agreement with ARPES 
results 



spectra in DMFT, see Fig. [31 Microscopically, infrared 
LFL behavior (narrow "Kondo" resonance in DMFT) is 
killed off by strong scattering between the Mott-localised 
and quasi-itinerant components of the (DMFT) matrix- 
spectral function, due to sizable U' , Jh, and is intimately 
tied with the OS insulator-metal transition in the five- 
band Hubbard model we use. Thus, our selective-Mott 
metal is a MO counterpart of the FL* metal 24 1 . 

If our proposal is to hold, a range of responses must 
find a consistent explication without additional assump- 
tions: we now show this is indeed the case, (i) A di- 
rect comparison between DMFT spectra and (AR)PES 
data show very good accord: in addition to describ- 
ing the overall PES lineshape very well, the DMFT 
spectra (inset of Fig. [3]) also resolve two peaks in the 
low-energy (—0.7 < U) < 0.0 eV) range seen in recent 
work for their "semiconducting" samples. Further, 
the renormalised FS is composed of predominantly xz, yz 
orbital states, and absence of hole pockets in LDA thus 
persists in DMFT, in accord with ARPES. Absence of 
LFL quasiparticles naturally implies broad ARPES line- 
shapes without LFL quasiparticles for electron-doped 
samples, as well as sizable SWT over an energy scale 
0(1.5 — 2.0) eV in polarised X-ray absorption (XAS) 
studies: this could be tested in future, (ii) Transport 
properties in DMFT are directly computable from the 
full DMFT propagators, since vertex corrections enter- 
ing the Bethe-Salpeter equation are small enough to be 
neglected in multi-band cases 25j. In Fig. |4l we show 
the dc resistivity for various n. Clear (Mott) insulating 
behavior obtains for n = 6.0, as expected. With 5 > 0, 
a doping-dependent crossover from a "high"-T insulator 
to a low-T (lower inset in Fig. [3} bad metal obtains. Ob- 



FIG. 4: Resistivities for insulating (solid curve) and electron- 
doped phases of Kj,Fe2-nSe2 as a function of T in good accord 
with transport data 0, H(| . 



viously, this crossover scale, marked by the maximum 
of pdc(T), increases with 5. The upper inset of Fig. [4] 
shows that small U yields metallic behavior up to high 
T, in conflict with data. At the very least, U > 3.5 eV 
is needed to obtain accord for 5 = 0, but correct insulat- 
ing behavior is only found for U = 4.0 eV. (iii) Further 
corroborative support comes from optical conductivity. 
In Fig. [5j we show the total optical conductivity, cr(w), 
within DMFT, using the multi-band version of the gen- 
eral DMFT formalism 25j. a(u>) rises with w at small u, 
and the optical gap, A ~ 0.3 eV, is in good accord with 
data [lo| . Further, as also seen, bad-metallic conduc- 
tivity, along with large-scale (> O(2.0) eV) transfer of 
dynamical spectral weight, obtains for 8 > 0. It is par- 
ticularly noteworthy that good accord is seen not only 
at low energy (0 < w < 0.3 eV), but up to 1.0 eV, as 
seen by direct comparison with data (inset to Fig. [5]). 
Comprehensive accord with PES and (dc and ac) con- 
ductivities thus provides strong support to the view that 
KFei.6Se2 is a Mott insulator (MI), and that the doping- 
driven transition to a bad metal must then be viewed in 
terms of a doped MI. 

These findings put strong constraints on mechanisms 
of SC. Since LFL quasiparticles are not stable excita- 
tions in the bad metal, instabilities to ordered states 
via (particle-particle) BCS-like pairing of well defined 
LFL quasiparticles are untenable. In non-LFL, bad met- 
als coherent one-electron hopping term is irrelevant in 
an RG sense [23}. Hence, ordered states can only arise 
via coherent two-particle hopping, which becomes more 
relevant when one-electron hopping term is irrelevant: 
the situation is analogous to coupled D = 1 Luttingcr 
liquids [2jJ. Given local incoherent metallic state(s) 
within DMFT, the above analogy tells us that resid- 
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FIG. 5: (Color online) Optical conductivity for insulating 
(solid curve) and electron-doped phases of K a: Fe2-j / Se2 as a 
function of energy. The good accord with experimental results 
(inset) up to Q ~ O(l.O) eV is seen. 



ual, inter-site and inter-orbital (in multi-band systems) 
two-particle interactions can generate ordered states di- 
rectly from the bad metal. This is the philosophy used 
earlier [28[ for the 1111-FePn systems. As in earlier 
work, restricting ourselves to the xz,yz orbital sector, 
the effective pair-hopping term (second order in t a b) 

is ~ -1/2 Efc,^,„,6 ^')4a/!,fc', ff ' C (>,f=^' C ^- 

Decoupling in the particle-hole and particle- 



particle channels gives H$ F = J2 k ,a,bl( A ab ' ( k ) c laa c kb* + 
h.c) + (A%(k)cl at cl kbi + h.c)}, where A$(fc) = 

{l(k)c J kbtT Ckaa) and A y ab ' (fc) = {^{k)c k baCka~ a ) with 
~f{k) =cosk x +cosk y + acosk x cosk y for the frustrated case 
of Fe-based systems. These represent orbital nematic 
(with orbital order and lattice distortion) [29[ and inter- 
orbital pairing 28] instabilities. Extending DMFT to 
study both these orders in K x Fe2- y Se2 is more prob- 
lematic, however: the large moment, [ip e = the 
block- spin moment M ~ and the block spin-AF 

order [30( suggests that both, possible orbital order [31 1 
with lattice distortion and subsequent AF, as well as 
SC instabilities must involve coupling between four-site 

(2) 

plaquettes, beyond what our H MF would give. Since 
the Mott transition already occurs at high T, a way to 
proceed might involve using the present DMFT results 
as a template for deriving an appropriate low-energy, 
plaquette-centered model using the active xz,yz orbital 
states to address these issues as recoginsed by Baskaran 
in a different approach [32j. This is currently underway, 
and will be reported in future. 

In conclusion, using LDA+DMFT for a minimally re- 
alistic five-band Hubbard model, we resolve the issue of 



a doped Mott- vs band insulator physics in K2;Fe2-ySe2 
systems in favor of the former, Mott view. Good quan- 
titative accord with key spectral and transport data 
in a sizably correlated picture confirms this view, and 
strongly suggests close underlying similarities (in spite of 
very different chemistry) between SC emerging here from 
a doped multi-orbital Mott-Kondo-insulating state with 
ci-wave SC in doped high-T c cuprates. 

L.C. thanks the Physical Chemistry department at 
Technical University Dresden for hospitality. 
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